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I 

Abstract 

The stochastic Euler scheme is known to converge to the exact 
solution of a stochastic differential equation with globally Lipschitz 
coefficients and even with coefficients which grow at most linearly. 
For super-linearly growing coefficients convergence in the strong and 
numerically weak sense remained an open question. In this article 
we prove for many stochastic differential equations with super-linearly 
growing coefficients that Euler's approximation does not converge nei- 
ther in the strong L p -sense nor in the numerically weak sense to the 
exact solution. Even worse, the difference of the exact solution and of 
the numerical approximation diverges to infinity in the strong L p -sense 
and in the numerically weak sense. 

1 Introduction 

An important numerical scheme for simulating stochastic differential equa- 
tions (SDEs) is Euler's method (see e.g. [T4l [20] ) (a.k.a. Euler-Maruyama 
scheme). If the coefficients of an SDE are globally Lipschitz-continuous, then 
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standard results (e.g. Chapter 10 and 14 in [T4]) show convergence of the 
Euler approximation in the strong and numerically weak sense to the exact 
solution of the SDE. It remained an open question whether the Euler ap- 
proximation also converges in the strong or numerically weak sense if the 
coefficients are not globally Lipschitz-continuous, see e.g. Section 1 in [8J. In 
this paper we answer this question to the negative. More precisely we prove 
for a large class of SDEs with super-linearly growing coefficient functions 
that both the distance in the strong L p -sense and the distance between the 
p-th absolute moments of the Euler approximation and of the exact solution 
of the SDE diverges to infinity for any p G [1, oo). Thus the Euler scheme 
does in the strong and numerically weak sense not produce an appropriate 
approximation of the exact solution of such an SDE. 

For clarity of exposition we concentrate in this section on the following 
prominent example (also considered in [HI [21], 23J). Let (X t ) t > be the unique 
strong solution of the one-dimensional SDE 

dX t = -X t 3 dt + o dW tl X = x , (1) 

where (Wt)t>o is an one-dimensional standard Brownian motion and a > 0, 
x G K. are constants. The Euler approximation (Y^)k=o,...,N of (Xt) te r 0T u 
T < oo fixed, is defined recursively through Yq* = xo and 

Y&i = Yu ~ | {Ykf + * (WW - Wq) (2) 

for every k = 0, . . . , N — 1 and every iV G N := {1, 2, . . .}. We will prove in 
Section [2] that the //-distance of the exact solution Xt and of the numerical 
approximation Y^ diverges to infinity as iV tends to infinity, that is, 

lim E\X T - Y$\ P = oo (3) 

for every p G [1, oo), see equation ffTTI) . Thus strong convergence of the Euler 
approximation fails to hold. In addition numerically weak convergence (see 
e.g. Section 9.4 in [H]) (not to be confused with stochastic weak convergence) 
fails to hold, that is, 



lim 

N—fOo 



E\X T \ P -E\Y$\ P 



oo (4) 



for every p G [l,oo), see equation ([TBI . Please note that this divergence of 
the Euler approximation is not a special property of equation ([TJ. We es- 
tablish this divergence for a large class of SDEs with super-linearly growing 
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coefficients in Section [3J Moreover our estimates are easily adapted to prove 
divergence of other numerical schemes such as the Milstein scheme. For sim- 
plicity we restrict ourselves to the Euler scheme. Presence of noise, however, 
is essential. In the deterministic case, e.g. (JTJ and (jSJ) with a = 0, the Euler 
scheme does converge and both ([3]) and (j3J) fail to hold. In addition note 
that - in contrast to strong and numerically weak convergence - pathwise 
convergence does hold if the coefficients are locally Lipschitz-continuous, see 
e.g. Gyongy (1998), Fleury (2005) or Jentzen et al (2009). 

Next we indicate why the L 1 -norm of Euler's approximation diverges if 
the growth of the coefficient functions is not linearly bounded. The Euler 
approximation of the deterministic equation, that is, (pQ) and (j2J) with a = 0, 
is known to be unstable if x is too large (see e.g. Chapter 6 in [3]). A central 
observation of this article is a quantification of this instability. For example, 
if Y N = x = N and T = 1, then 

Yf — N — —N 3 = N-N 2 « -N 2 (5) 

and therefore 

Y 2 N = Y* - ^ (Y^) 3 « (-N 2 ) - ^ (-N 2 ) 3 = -N 2 + N 5 ^N 5 >N A 
and 

Y 3 N = Y 2 N -^ {Y 2 N f « (iV 5 ) - L (ATS) 3 = iys _ AT14 w .^14 < _at8 

for every large iVGN. Iterating this calculation, yields l^j^l > iv( 2fc ) for ev- 

ery k — 0, 1, . . . , AT and, in particular, \Y$\ > n( 2N ) for every large iVGN. 
Thus grows double-exponentially fast in N. Now in the presence of noise 
(er > 0) there is an exponentially small event that the Brownian motion 
leaves the interval [— 2A^, 2N\. On this event the deterministic dynamics lets 
the approximation grow double-exponentially fast. Consequently as being 
double-exponentially large over-compensates that the event has an ex- 
ponentially small probability, the L 1 -norm of the Euler approximation is 
unbounded in N. 

Having sketched why the Euler scheme diverges for many SDEs whose 
coefficient functions are not linearly bounded we take a short and incom- 
plete look into the literature where we concentrate on the non-globally Lip- 
schitz case and on relatives of the Euler scheme. It is a classical result 
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(e.g. [HI [2D]) that the Euler approximation converges in case of globally 
Lipschitz-continuous functions. In contrast to our results on super-linearly 
growing coefficients, Yan (2002) proves numerically weak convergence of the 
Euler scheme if both drift and diffusion function have at most linear growth. 
Higham et al (2002, Theorem 2.2) provide an instructive condition for strong 
convergence saying that local Lipschitz-continuity of the coefficients together 
with boundedness of a p-th moment, p > 2, of the approximation and of the 
exact solution suffices for strong L 2 -convergence. Of course our results sig- 
nify that the p-th moment of the explicit Euler approximation is unbounded 
if the coefficients grow super-linearly and thus this condition does not apply 
in that case. Zhang (2006) and Berkaoui et al (2008) prove strong conver- 
gence for a class of drift and diffusion functions with a singularity. Yuan and 
Mao (2008) obtain the rate of convergence of the stochastic Euler scheme for 
locally Lipschitz-continuous coefficients if these grow at most linearly. An ex- 
plicit bound for the strong L 1 error is established in Higham and Mao (2005) 
for the mean-reverting square-root process (a.k.a. CIR process), which has 
linearly bounded coefficients. Bernard and Fleury (2001) establish conver- 
gence in probability under weak hypotheses like local Lipschitz-continuity. 
Pathwise convergence results can be found in Gyongy (1998), Fleury (2005), 
Jentzen and Kloeden (2009) or in Jentzen et al (2009). 

A number of authors obtain convergence for modified Euler schemes. 
Milstein and Tretyakov (2005) consider a modified Euler scheme for non- 
globally Lipschitz coefficients. They obtain numerically weak convergence 
by discarding trajectories of the Brownian motion which leave a sufficiently 
large sphere. As we have seen above, these trajectories are responsible for 
the divergence of the absolute moments. Lamba et al (2007) prove strong 
convergence of an Euler scheme with an adaptive time-stepping algorithm 
for locally Lipschitz-continuous coefficients. A completely different adaptive 
time-stepping algorithm with focus on long time approximation is proposed 
by Lemaire (2007). For results on implicit Euler schemes consult e.g. Hu 
(1996) or also Talay (2002) and the references therein. Finally, Higham et 
al (2002) showed strong convergence of an implicit Euler scheme via the 
instructive condition cited above. 

Simulations do not always reveal the divergence of the Euler approxima- 
tion. See Table [1] for simulations in which the approximate value is within 
distance two of the exact value. To understand this consider the determin- 
istic case, that is, ([1]) and (j2J) with a — 0. Then the approximation explodes 
if \xq\ > ^2N/T. Consequently in the case a = 1 = T it is likely that 
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the approximation diverges if the Brownian path reaches level 
probability hereof is bounded above by 



/ 2N. The 



P 



sup \W t \ > V2N 

0<i<l 



< 4-P 



< 



e~ N < 2e~ N 



(6) 



for every N G N, where we have used the reflection principle and a standard 
estimate, see e.g. Theorem 21.19 in [13]. Thus the probability that we see 
an explosion of the approximation is negligible for N > 100, say. Even more 
if the simulation is repeated N 2 times, say, then the probability to see an 
explosion is bounded above by about 2N 2 e~ N for every N G N, which is again 
negligible for N > 100, say. Therefore even in a large number of simulations 
it is unlikely to encounter an approximation which diverges. 

The rest of this article is organized as follows. To clearly present the 
idea of this article, we exemplify the divergence of the Euler scheme for the 
SDE (pQ) with a = 1, x = in Section [2] (Theorem [Q. A general result 
(Theorem [2]) and several examples for the divergence of the Euler scheme 
are provided in Section [3l Simulations in Section H] illustrate this divergence. 
The proofs are postponed to the final section. 



2 A simple counterexample 

Fix T > and let (Q, J 7 , P) be a probability space on which a scalar standard 
Brownian motion W : [0, T] x Q — > R with continuous sample paths is defined. 
Then, in this section, we consider the SDE (pQ) with a — 1 and Xq = 0, i.e. 

dX t = -X? dt + dW t , X = (7) 

for every t G [0,T]. More precisely, let X: [0,T] x Q — > R be the pathwise 
unique stochastic process with continuous sample paths, that satisfies 

X t (uj) = - [ (X s (u)f ds + W t (u) (8) 

for every t G [0, T] and every w G fi (see e.g. Theorem 2.4.1 in [TB] for 
existence and uniqueness of such a process). Then the Euler scheme denoted 
by jyi3(M)-measurable mappings : -> R, k = 0, 1, . . . , N, N G N, is in 
the setting of equation (jTj) given by Y^(u) = and 

Y&M = Y k N H - 1 (^(a;)) 3 + (mW(u,) - W,_r(u)) (9) 
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for every k = 0, 1, . . . , N — 1, every N G N and every u G Vt. We now show 
the divergence of Euler's method for the SDE (J7j). 



Theorem 1. Let T > 0, the solution X of (151) and £/ie Euler approximation 

-A 



y.tf, Jfe = 0, 1, . . . , iV, iV G N, 6e as a&ove. Taen 



lim E |X T - Y£ = oo. 



AT— »oo 



(10) 



Remark: Due to Jensen's inequality Theorem [T] in particularly implies 

lim E [\X T - Y$\ p ] = oo (11) 



N^oo 



for every p G [1, oo). In addition since the exact solution X : [0, T] x fi — > R 
satisfies E |X T | P < oo for every p G [1, oo) (see e.g. Theorem 2.4.1 in [18J) we 
get 



(E|1^| P )^ = (E\Y»-x T + x T \y 

> (e\y£-x t \ p )* - {E\X T \ p )p 



(12) 



oo 



as N — ► oo due to equation IjTTj) . Therefore we obtain 



lim |E|X t | p -E|Y^ |P| 



OO 



(13) 



for every p G [1, oo). 



Proof of Theorem\J\ First of all we define r^ := max (3-^, l) for every N G 
N. Then we consider the sets 



4 JV 



| w G ft 


sup 




fe=l,...,JV-l 



W{*+i)r (a;) — Wkr (ui) 



< 1, 



for every iVeN. Note that Q N G J 7 for every iV G N. 
Let now N G N be arbitrary. Then we claim 



Y k N (u)\>(r N ) 



(2< fc - 1 )) 



V wen 



iV 



(14) 
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for every k = 1, 2, . . . , N. We prove equation (1141) by induction on k = 
1, 2, . . . , N. In the base case k = 1 we have 



Y»(u>) 



Wt M - W a (u) 



(15) 



for every uj G Qn due to the definition of the Euler scheme (see equation 
and due to the definition of the set Qn- Therefore we assume that equation 
(114"]) holds for one k G {1, 2, . . . , iV — 1}. In particular we obtain 



ln%oi>M 2 >r„>i 

for every u G fijy Due to the definition of the set Ojv, we also have 



W(k+i)T (u>) — WkrUjj) 



< 1 



for every uj G Qn- Hence, equation (jHJ) and equation (1T?1) imply 



> 



I «V)) 

T 



- {Y k N (uj)) 3 + (Ws^t{u) - W ¥ (oo) 

- \ Y k N (uj)\ - W(k + i)T (u) - Wkr(u) 



for every u; G fi/v- Therefore equation ( fTBl) implies 



TV 



T 



2 in^H 



'AT/ 



N 

T 



>\Y k "(u;)n-r N -2)>\Y»(u 



(16) 



(17) 



(18) 



(19) 



for every oj G fijv, since ^-r^r — 2 > 1 by definition. Therefore the induction 



TV 

hypothesis, i.e. that equation (fT4"]) holds for k, implies 

2 /„ 



(20) 
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for all oj G Qn, which indeed shows that equation (1141) holds for k + 1. Hence, 
we have shown equation (j!4p for every k — 1,2, ... ,N and since N G N was 
arbitrary, we obtain 

\Y»(u)\>(r N )( 2(N - 1} ) (21) 
for every u G Qn and every iVeN. This yields 



E|y*n 



\Y$(u)\F(du) > I \Y$(u)\P(dtv) >F[Q N ] ■ 2( 2<A " 1) ) (22) 
for every iVeN. Since 
F[Q N ] 



p 


sup 


W(k + 1)T — WkT 


< 1 


•P 




Wt 


> r N 




_jfe=l,...,JV— 1 


N N 








N 





> p 



p 



> p 



> 



1 



sup \W t \ < - 

0<t<T 4 

sup \W t \ < r 

0<t<T * 

sup \W t \ < - 

0<t<T 2 

p 



■ p 

• p 
1 

' 4 



Wt 

N 

N 
T 



> r N 



Wt 

N 



(23) 



sup \W t \ < 



0<t<T 



■ e t 



N (r N f 



for every AfeN due to Lemma [T] in Section below, we obtain 



1 



•P 



4VT 

for every N G N. This shows 



sup \W t \ < - 



0<t<T 



_N_ 

e T 



(24) 



lim ElY^I > 



4VT 



•P 



sup \Wt\ < 

0<t<T 



1 




< - 


• lim ( 


1 - 2 





4VT 



•P 



sup \W t \ < 



0<t<T 



• lim ( e V^J ■ 2 

N^oo 



(r N f .2( 2(Ar - 1) ) , 

' I = OO. 



Since, as mentioned above, the exact solution X : [0, T] x f2 — > R of the SDE 
(J7j) satisfies E \Xt\ < oo (see e.g. Theorem 2.4.1 in [15]). we obtain 



lim E |X T - Y^l > lim E \Y£ I - E |X T | = oo 



N I 



jV-»oo 

which is the assertion. 



N^oo 



(25) 
□ 
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3 Further counterexamples 

Throughout this section assume the following setting. Fix T > and let 
(Q, T, P) be a probability space with normal filtration (J~t) te r r i • In addition 
let W : [0, T] x Q — > R be a scalar standard Brownian motion with respect 
to {^Ft)te[o,T\ w ith continuous sample paths and let £ : fi — > R be a JF /£>(R)- 
measurable mapping. Let f,g: R — > R be two i3(R)/i3(R)-measurable func- 
tions such that the SDE 



dX t = f(X t )dt + g(X t )dW t , Xo = £, t€[0,T] 



(26) 



has a solution X : [0, T] x f2 — > R. More precisely we assume that X : [0, T] x 
Q — > R is a predictable stochastic process satisfying 



P 



/ T |/(X,)| + |s(X s )| 2 d S < 



and 



P 



X t = i+ f f(X s )ds+ [ g{X s 
Jo Jo 



oo 



1 



(27) 



(28) 



for every t G [0, T]. Then Euler's approximation for the SDE fl2"B"l) - denoted 
by J^/S(R)-measurable mappings Y^: fi -»• R, k = 0, 1, . . . , N, N G N - is 
given by Yq^u;) = and 

= + {Y k N (u)) (Ws^r( U ) - W,*{uj) (29) 

for every k — 0,1, ... , N—l, every iVeN and every G fi. Now we formalise 
the main result of this article which asserts the divergence of Euler's method 
for the SDE f[2"6l) with super- linearly growing coefficients. 

Theorem 2. Assume that the setting above is fulfilled with P [#(£) 7^ 0] > 
and letC>l,f3>a>l be constants such that 



max(|/(x)| , \g(x)\) > 



\x\ 



C 



min (1/(^)1 , \g( x )\) < c \ x \ a 



(30) 



for all \x\ > C . If the exact solution of (126]) satisfies E \Xt\ v < 00 for one 
p G [1, 00), then 



lim E \X T — Y, 



N\P 

N 



00 and lim E|Xrl P — \Y M 

N^oo 1 1 1 1 JV 



00, 



(31) 



where Y k , k = 0, 1, . . . , N, N G N, is defined in equation f[2"9"j) . 
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Condition (1301) should be read as follows. Either the drift function or the 
diffusion function grows in a higher order than linearly and the remaining 
function grows slower than that. More formally it suffices to show that either 



l/(*)|>^-, \9(x)\<C\x\ a (32) 

for all \x\ > C or 

\g{x)\>^, \f(x)\<C\x\ a (33) 

for all \x\ > C and for some constants (3 > 1, (3 > a > 0, C > 0. Please 
note that our estimates need (3 > 1. A super-linear drift function such as 
f(x) = xlogx is too small for our estimates. The assumption that the 
diffusion function does not vanish on the starting point ensures presence of 
noise in the first time step. 

The proof of Theorem [2] is deferred to Section [5j In the remainder of 
this section we apply Theorem [2] to a selection of examples. Please note that 
the coefficients in the following examples are locally Lipschitz-continuous 
and satisfy an appropriate one-sided Lipschitz condition. Therefore the p-th 
absolute moment in each example is finite for every p G [1, oo) according to 
Theorem 2.4.1 in [18]. For every example we check that assumption (!30|) is 
satisfied. Thus both the distance in the strong L p -sense and the distance 
between the p-th absolute moments of the Euler approximation and of the 
exact solution diverges to infinity for any p 6 [1, oo) in each of the following 
examples. 

The introductory example The example in Section [1] is 

dX t = -Xfdt + adW t , X = x E R (34) 

where a > 0. The dominating coefficient is the drift function with dominating 
exponent (3 = 3. The diffusion function has exponent zero and we may choose 
a = 0. The constant C in equation fl32l can be chosen as C = max(l, a). 

Stochastic Ginzburg-Landau equation The Ginzburg-Landau equa- 
tion is from the theory of superconductivity. It has been introduced by 
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Ginzburg and Landau (1950) [5] to describe a phase transition. Its stochas- 
tic version with multiplicative noise can be written as 

dX t = ^(77 + ^o- 2 )X t - XX^j dt + aX t dW t , X = x E (0, 00) (35) 

where rj > 0, A, a > 0. Its solution is known explicitly (see e.g. Section 4.4 
in p3!) 

Xt = x exp (rjt + aWt) ^ 



'l + 2x 2 X f* exp (2r)S + 2aW s ) ds 

for t > 0. Since 



a 



7] + — J X — Xx 



> A |x| 3 - \t) + —J \x\ 

A I 1 3 A , .3 / (J 2 

A , ,3 . A / 2 2ri + a< 
\x\ H — \x\ \ \x\ — 



a 



(37) 



2 1 1 2 1 1 V A 

> — X H X X : > — X 

2 1 1 2 1 1 V A / _ 2 



for all |x| > max ^1, 2,? ^ cr J , the constants in equation (|32|) can be chosen 
/? = 3, a — 1, C = max ^1, a, |, 



as 



2 2>7+o- 2 
A 

Stochastic Verhulst equation The Verhulst equation is an ordinary dif- 
ferential equation and is a simple model for a population with competition 
between individuals. Its stochastic version with multiplicative noise can be 
written as 

dX t = ((rj + \° 2 )X t - XX^j dt + aX t dW t , X = x E (0, 00) (38) 

where rj, X, a > 0. Its solution is known explicitly (see e.g. Section 4.4 in p~4] ) 
and is given by 

x t = *o*v(yt + °Wi) (39) 

1 + x A f exp (r/s + aW s ) ds 

for t > 0. The dominant exponent of the drift function is two and of the 
diffusion function is one. Thus we may choose (3 = 2 and a — 1. The 

constant C in equation (1321) can be chosen as C = max (a, |, 2ri+ ^ J . 
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Feller diffusion with logistic growth The branching process with logis- 
tic growth (see e.g. [IB]) is a stochastic Verhulst equation with Feller noise. 
It solves 

dX t = \X t (K-X t )dt + (Ty/xldW t , X = x e(0,oo) (40) 

where \,K,a > 0. There is no explicit solution for this equation. However 
it features the following self-duality 

E*exp (~X t y) = E^exp ( ~xX^\ , Vt,x,y> 0, (41) 



x a 2 J \ a 2 7 

where K x refers to the starting point Xo = x, see [10J. As constants for 
equation (1321 serve f3 — 2, a — 1, C — max(l, |, a, 2K). 

Ohta-Kimura model The Ohta-Kimura model ([22]) is a population mo- 
del with two types, fixed total population and random selection. Its diffusion 
equation is 

dX t = aX t {\ - X t ) dW t , X = x E (0, 1) (42) 
where a > 0. Since 



\ax(l — x)\ > a \x\ —a\x\ = — \x\ H — \x\ —a\x\ 
iv 1 1 2 1 1 2 1 1 11 

° I |2 , ° I I l\ I w ff I |2 

= — \x\ H — \x\ { \x — 2) > — x 
2 i i 2 i i vi i 2 



(43) 



for all \x\ > 2, Theorem [2] applies with (3 = 2, a = and C = max(2, -) in 
equation (153^) . Note that this is an example with the diffusion function being 
the dominant coefficient. 



4 Simulations 

In this section we present simulations which illustrate the dichotomy between 
the double-exponential growth of the deterministic dynamics (for some initial 
values) and the exponentially small probability of the event that the noise 
is sufficiently large. For this purpose we choose an equation with an explicit 
solution to compare with. Consider the Stratonovich equation 

dX t = -X? dt + aX t o dW u Xo = 1, t e [0, 3]. (44) 
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This equation agrees with the stochastic Ginzburg-Landau equation (1331) in 
the case rj = 0, A = 1, x = 1. Its explicit solution is given by equation (1551) 
with 77 = 0, A = 1, Xo — 1. 

Table H] shows Monte Carlo simulations of the second moment both of 
the explicit solution and of the Euler approximation for five values of a. 
For each value of a the table contains ten simulation runs. The number of 
time steps is N = 10 3 . In the case a = 7 the parameters of the simulation 





a = 2 


a = 4 


(7 = 5 


(7 = 6 


a = 7 


Exact solution E (X3) 2 


0.4739 


0.9598 


1.2357 


1.5423 


1.8861 


Simulation 1 


0.4577 


0.7448 


0.7243 


NaN 


NaN 


Simulation 2 


0.4549 


0.7518 


0.7094 


0.5097 


NaN 


Simulation 3 


0.4632 


0.7330 


0.6912 


NaN 


NaN 


Simulation 4 


0.4580 


0.7544 


0.7150 


0.5378 


NaN 


Simulation 5 


0.4619 


0.7458 


0.6998 


0.5197 


NaN 


Simulation 6 


0.4610 


0.7664 


0.7327 


0.5243 


NaN 


Simulation 7 


0.4606 


0.7542 


0.7189 


NaN 


NaN 


Simulation 8 


0.4690 


0.7451 


0.7122 


NaN 


NaN 


Simulation 9 


0.4667 


0.7336 


0.7146 


0.5475 


NaN 


Simulation 10 


0.4591 


0.7605 


0.6659 


NaN 


NaN 



Table 1: Simulation for the SDE (H}. The number of time steps is N = 10 3 . The 
number of Monte Carlo runs for the approximation is 10 5 and for the exact solution 10 7 . 



are such that the Euler approximation produces the value "NaN" (NaN is 
the IEEE arithmetic representation for "not-a-number" , here this is due to 
an operation "Inf — Inf" where Inf is the IEEE arithmetic representation 
for positive infinity). In contrast to this the Monte Carlo simulations in 
the cases a e {2,4,5} are all finite. As explained in the introduction the 
Euler approximation grows double-exponentially fast only on an event of an 
exponentially small probability. For a G {2,4, 5} this probability is so small 
that none of the simulation runs produced a large trajectory. The last but 
one column in Table [T] exemplifies a parameter setting in which the event of 
large growth has a probability such that in some Monte Carlo simulations no 
explosion occurs and in some Monte Carlo simulations at least one explosion 
in 10 5 runs occurs. 

The values in Table [T] are either within radius two of the true value or 
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"NaN". This is due to the double-exponential growth of the deterministic 
system for some initial values. If the simulation starts to grow, then it 
reaches "NaN" very quickly. We encountered similar behaviour for other 
exponents greater than two. In order to see double-exponential growth of the 
Euler approximation in a plot we consider an exponent close to one. More 
precisely we plot the Monte Carlo simulation of the first absolute moment of 
the Euler approximation of X w where 



Absolute moment E|Y^J| of the Euler scheme versus N=1,2,3,...,50 




5 10 15 20 25 30 35 40 45 50 
Number of time steps N 



Figure 1: First absolute moment E \ Y$\ of the Euler scheme for number of time steps 
AT = 1,2, 3, ... ,50. Number of Monte Carlo runs: 10.000 

dX t = -lOsgnpQ) |X t | la fit + AdW t , X = 0, (45) 

see Figure [1] (see also Table [2]). Please note that the graph resembles an 
exponential function and that the y-axis is logarithmic. Thus the growth of 
the graph in Figure [T] is indeed close to a double-exponential function. In 
addition we should mention that some fine-tuning was necessary to obtain 



14 



N 




N = 1 


10.1160 


N = 2 


444.5797 


N = 3 


1.3649 • 10 4 


N = 4 


3.7566 • 10 5 


N = 5 


1.0295- 10 7 







N 


E ^ 


N = 48 


3.4449 ■ 10 260 


iV = 49 


2.6631 ■ 10 271 


iV = 50 


6.1140- 10 297 


N = 51 


7.2342 ■ 10 3U1 


A^ = 52 


Inf 


N = 53 


Inf 



Table 2: Simulation values of the first absolute moment E \ Y^\ for N = 1, 2, 3, ... , 53 

suitable parameters for which the simulated absolute moment grows but is 
not "NaN". 



5 Proofs 

Lemma 1. Let (Q, J 7 , P) be a probability space and let Z : Q —>■ K be an 
T /B(R) -measurable mapping which is standard normal distributed. Then we 
have 



P 



\Z\ > x 



> 



x e~ 



P 



\Z\ e [x,2x] 



> 



xe 



-2x z 



for every x G [0, oo). 

Proof of Lemma El We have 



(46) 



P[|Z| > x] = 2 ■ F[Z > x 



2 I / -^=e~^ds 
L V2n 



> 2 



1 -*-A 

e 2 as 



> 2 (V2- 



l _(y^r 

— e 2 



2tt 



'2(V2-1)\ _ x 2^.n 
1 xe > 



'2ti 



(47) 
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for every x G [0, oo). Moreover, we have 

F[\Z\e[x,2x]] = 2-F[Z e [x,2x]] = 2 

1 



2;r 



-e 2 ds 



2tt 



> 2x 



^7T 



(Ml 
=e ' 2 



>2tx 



2 x e 



xe- 2x ~ > 



for every a; G [0, oo) 



(48) 



□ 



Proof of Theorem^ First of all we have P[ \g{£,)\ > 0] > 0. Therefore there 
exists a real number K > 1, such that 



/i := P 



is(oi iei+T|/(oi<K 



> 0. 



Then we define 



(49) 



/'a :-iuax(2,C, ( — + 2C 2 



G [C, oo) 



(50) 



and consider the sets fl^ G T given by 



Q N := < uj E Q 



W(k+i)T (lu) — Wkr (a;) 





"T 


2T" 


G 


_iV' 





VJfe = 1,...,JV- 1, 



WMw) - W (u) 

N 



>K(r N + K) 



I^H)I>^, Uu)\+T\f{t-{u))\<K\ (51) 



for every iV G N. Note that 



T {rtl) ^>^(^E + 2C^>2 + 2TC 



NC 



NC \ T 



N 



(52) 



for every N G N due to the definition of r N , N eN. 
Let now iVGNbe arbitrary. Then we claim 



|^)|>M( a(M ) 



V weO 



(53) 
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for every k = 1, 2, . . . , N. We prove equation f)53p by induction on k G 
{1, 2, ... , N}. In the base case fc = lwe have 



= ch + -/ (e(w)) + (e(w)) (wt ( W ) - w (w)) 

>\g(m)\ 
>\g(m)\ 



Wt(uj)- W a (u) 



\m\-^\f(m)\ 



Wt m - W n (u) 



(54) 



K 



> 



K 



Wt(u)-W (u) 



- K > K = r N 

A 



for every uj G Qn due to equation ( 129|) and equation ( loTT) . For the in- 
duction step -> fc + 1 we assume that equation ( 153|) holds for one G 
{1, 2, . . . , N — 1}. In particular we obtain 



\Y k N {u)\>r N >C 
for every u E Qn, since rjy > 2. We also have 



(55) 



T 



> 



T 



iV 
> max 



f(Y h N (u)) + g(Y h N (u)){\\' l , tlll U<} t%M 
T 



mm 



N 
T 
N 



f(Y»(uj)) , \g(Y k "(u))\ ■ W w (u) - W*r{u) 



f{Y»{u))\Ag{Y»{Lo)) ■ W (k+1)T (u)-W*r(u) 



and 



|n+l(^)| > max (L |/(lfH)| , § 



T 



2T 



iV 



iV 



T 



(56) 



>-max(|/(F^M)|,|^M)|) 
2T 

- — min {\f(Y k N (u))\ , |^H)|) - |n 



OJ 
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for every uj G Vt^ due to equation (1251) and equation (1ST]) . Therefore equation 
(157)]) and equation (150]) yield 



In+iMl > 


T 

ATC ' 


Y k N i 




/3 


2TC | 

A" 






a j 


Y k N H\ 


> 


T , 

A^C 1 


Y k N ( 


»| 




2TC - 

AT 1 


Y k N i 


H 


a j 


Y k N {u)\ a 




|^( 


»)|° 


/T 
V A^C 


\Y k N (u> 


)t 


-a) _ 


2TC \ 
N ) 



(57) 



> m^r m ( ^ q) - ^ - 1) > m»r 

for every u E where the last step follows from equation (l52j) . Hence the 
induction hypothesis yields 

IO-)l > |n»r > (W^ 15 ))" = M (flfc) (58) 

for every ui G Ojv, which proves equation (153]) for + 1. Therefore equation 
(155]) indeed holds for every G {1, 2, ... , AT}. 

In particular - since iVeN was arbitrary - we obtain 

I^MI^M^ 1 ') (59) 
for every to G fijv and every iVGN. Thus 

E|r^| = / \Yg(u)\ F(du) > j \Y*{uj)\ F(dou) 

> / (r N )( a{N ^) ndco) = P M • M^- 1 ') > P ■ 2(* (iV - 1, ) 
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for every N G N. Hence we obtain 

W(k+1)T - W 



e\y£\ > p 



p 



fcT 
N 





T 


2T~ 


e 







Vk = l,...,N-l 



Wt - W 



>K(r N + K)] ■ 11 ■ 2( q(N 15 ) 











"T 


2T~ 




\ (iV-l) 












e 








] -P 




Wt 






JV 




_iV' 


~W_ 








N 



>K{r N + K) 



■ / ,.2K iV - 1) ) 
> | P 



(60) 




Wt 

N 





2VT 




) 











N 



P 




Wt 

N 



IN 



> y —K (r N + K) 



. / ,. 2 H iV - 1) ) 



for every iVeN, where // > is defined in (|49|) . Hence, Lemma [T] yields 



2 V iV 

9 (a< Ar - 1 ) > 

' 4 



A? 



> 



1 /t" 

2V 7v 1 ( 



4^ 



(61) 



> (L/I) . e -f *w*)\ ^ . 2 ( 

- I 2 V AW ^/T 



A? 



2T 



for every IVeN. This shows 

limE|r^|> lim fpffijvl ■2( a(JV_1) ) N ) 

AT-+oc 1 1 AT->oo V / 



/ie- 2T 
> • lim 

4\/T at-*oc 



OO. 



Due to Jensen's inequality we conclude 



lim E \ Y, 



N^oo 



N \P 

N 



OO. 



(62) 
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Hence, since E \X T \ P < oo, we have 

lim |E|X T | P -E|y^n = oo, lim E \X T - Yfi\ p = oo, (63) 

N—>oo N— >oo 

which is the assertion. □ 
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